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ABSTRACT 

We study thermal instability in a magnetized and partially ionized plasma with 
charged dust particles. Our linear analysis shows that the growth rate of the unstable 
modes in the presence of dust particles strongly depends on the ratio of the cooling 
rate and the modified dust-cyclotron frequency. If the cooling rate is less than the mod- 
ified dust-cyclotron frequency, then growth rate of the condensation modes does not 
modify due to the existence of the charged dust particles. But when the cooling rate 
is greater than (or comparable to) the modified dust-cyclotron frequency, the growth 
rate of unstable modes increases because of the dust particles. Also, wavenumber of 
the perturbations corresponding to the maximum growth rate shifts to the smaller 
values (larger wavelengths) as the cooling rate becomes larger than the modified dust- 
cyclotron frequency. We show that growth rate of the condensation modes increases 
with the electrical charge of the dust particles. 
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1 INTRODUCTION 

One of the most important dynamical processes in astro- 
physical plasmas is thermal instability. A detailed analysis 
of thermal instability in the linear regime was given in a 
well-know paper by Field (1965), and in large number of 
subsequent works (e.g., Heyvaerts 1974; Balbus 1986; Bal- 
bus & Soker 1989; Ibanez & Escalona 1993; Steele & Ibanez 
1999; Burkert & Lin 2000; Dib, Burkert & Hujeirat 2004; 
Hennebelle & Audit 2007). The form of net cooling func- 
tion is determined by the specific mechanisms of heating 
and cooling and it also depends on the degree of ionization 
of the gas and so, a more complete analysis of the thermal 
instability would require a corresponding equation for the 
degree of ionization (Defouw 1970; Goldsmith 1970). 

Thermal instability of a plasma in a uniform magnetic 
field has also been studied by Field (1965). It is obvious that 
for perturbations with wave vectors parallel to the field, the 
field does not affect the dynamics. If the wave vector is nor- 
mal to the magnetic field, the magnetic field reduces the 
growth rate of the unstable modes. Thus, a purely trans- 
verse magnetic field can prevent the thermal condensation 
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(e.g., Oran, Mariska & Boris 1982; Loewenstein 1990; Pio- 
ntek & Ostriker 2004). Hennebelle & Perault (2000) studied 
magneto-thermal instability analytically and numerically in 
the general case corresponding to the neutral interstellar 
medium. David & Bregman (1989) studied the effects of heat 
conduction and magnetic fields on the growth rate of ther- 
mal instabilities in cooling flows. Recently, Fukue & Kamaya 
(2007) revisited the effect of the ion-neutral friction of the 
two fluid on the growth of the thermal instability. 

It is also well-known that dust particles constitute an 
ubiquitous and important component of many astrophysi- 
cal plasmas including interstellar medium, stellar and plan- 
etary atmospheres, planetary nebulae and giant HII regions 
(e.g., Draine 2003). Dust particles can alter dynamical and 
thermodynamical properties of the plasma by transform- 
ing part of thermal energy into radiation (e.g., Kopp & 
Shchekinov 2007; Bora 2004; Ibanez & Shchekinov 2002). 
Birk & Wiechen (2001) have considered the radiative con- 
densation modes in a dense dusty plasma, where electrons 
are completely replaced by massive dust particles. However, 
in their analysis dust particles were considered only from the 
point of view of electrostatic interactions, while cooling pro- 
cesses were treated as independent of the presence of dust. 
Shukla & Sandberg (2003) studied radiation-condensation 
instability in a self-gravitating dusty astrophysical plasma. 
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Kopp & Shchekinov (2007) showed that positively charged 
dust particles strongly destabilize perturbations. 

In our model, neutrals are the thermodynamically ac- 
tive component and providing conditions for radiative cool- 
ing of the system, while ions, electrons and dust particles 
are the passive components. We show that when charged 
dust particles do not contribute to the cooling of the sys- 
tem, thermal stability properties of the system can modify. 
In the next section, we present general formulation of a mul- 
tifluid system including charged dust particles. Analysis of 
the unstable modes are given in section 3. We conclude by 
a summary of the results in the final section. 



2 GENERAL FORMULATION 
2.1 Basic multifluid equations 

For the sake of simplicity we will describe in what follows the 
dust component as an ensemble of particles of equal masses, 
ma, and equal charges, i.e. Ze. Our basic equations and 
the main assumptions are similar to Pandey & Vladimirov 
(2007). We take account of the different bulk velocities and 
densities of the neutral, electrons, ions and charged dust 
particles in the system. The continuity equation is 



^ + V.(*v,) 



0. 



(1) 



where pj and Vj are the velocity of the various plasma com- 
ponents and the neutrals, respectively. 
The momentum equations are 



,_, Vj x B . 

= -q 3 n 3 (E H ) - pjVjnVj, 

c 

pn ^~w~ + Vn - Vv ") = -V-P+ y^pj^wVj- 

e,i,d 



(2) 



(3) 



Note that velocities Vj are written in the neutral frame and 
E' = E + v n x B/c is the electric field in the neutral frame. 
j stands for electrons (q e = — e), ions (g; — e) and dust 
(qd = Ze), where Z is the number of charge on the grain. 
The other physical variables have their usual meanings. 

Also, the collision frequencies is (Draine et al. 1983; 
Draine 2003) 

< av > jn 

fin — ; Pn, 

rrij + m n 



(4) 



where < av >j„ is the rate coefficient for the momentum 
transfer by the collision of the j th particle with the neutrals: 

< av > in = 1.9 x 10" 9 cm 3 s"\ (5) 
T 



± I 

< av > en = 4.5 x 10" 9 ( — )5 C m 3 s" 1 , 



(6) 



and for small grains, we have < av >dn~< av > in , but for 
grains ranging between a few Angstrom to a few microns 
(Nakano & Umebayashi 1986) 



<av> dn =2.8x 10 5 (^)^ 



10- 5 cm' 



(7) 



where a is the grain radius. In our calculation, for the ion 
mass and mean neutral mass we adopt m, = 30m p and 



1.67 x 10~ 24 g is the proton 
mass. 

Defining the mass density of the bulk fluid and the bulk 
velocity as p w p n and u « v„, equations (1), (2) and (3) 
give the continuity and the momentum equations for the 
bulk fluid as 

g+V.(pu)=0, (8) 

p(^ + u.Vu) = -VP+i^. (9) 

The next simplifying assumption is that electrons and 
ions are assumed well coupled to the magnetic field which 
implies (3 e 3> /?» >• 1, where /3j = uj C j/uj n is the ratio of 
cyclotron u> c j = qjB/m,jC to the collision frequencies. Based 
on this assumption and using quasi-neutrality condition, wc 
have (Pandey & Vladimirov 2007; Ciolek & Mouschovias 
1993) 



where 

e = [i + ^]/3i 



(11) 



Note that equation (11) is valid only for the radial compo- 
nent. General expression for 6 contains Hall and Pedersen 
conductivities (e.g., Ciolek & Mouschovias 1993). Thus, the 
induction equation can be written as (Pandey & Vladimirov 
2007) 



= V x [(u x B) - ^-^-J x B] 

dt LV ; Zen d 1 



(12) 



We assume that dust particles do not contribute in cool- 
ing of the system. One should stress, however, that at tem- 
perature 3 x (10 7 — 10 8 )K dust particles are efficiently de- 
stroyed in collisions with the ions and electrons (Draine & 
Salpeter 1979), and therefore only in the temperature range 
between 10 6 K and 3 x 10 7 K does dust contribute sufficiently 
to the net cooling. So, we write energy equation as 

i-^i?i + ^- v -^ v ' T+ ^ v ^ - °' (13) 

where the total time derivative is d/dt = d/dt + u.V, and 7 
is the ratio of specific heats. Also, £1 represents the energy 
losses minus energy gains per unit mass. The coefficient of 
thermal conductivity K has the values K\\ and K± in di- 
rections parallel to and perpendicular to the magnetic field 
B. 

The net cooling function Q is written depending on the 
physical conditions of the systems. For example, we are in- 
terested in the structure formation in a typical H I dusty 
region with C II cooling. We assume that To = 100 K and 
K = 4.9 x 10 3 ergs cm -1 s _1 KT 1 (Ulmschneider 1970). 
Then, C II cooling is dominant when the temperature is 100 
K, according to Wolfire et al (1995). Thus, 



n = 2.54 x 10" 14 A c /cii(7 n H o +j e n e ) 

77lH 

x exp(-92/T)ergs -1 g~\ 



(14) 



where the fluid consists of H atoms with mass mu- Also, nx 
is the density of element X, the constant /en is the fraction 
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of C I in C II, and A c = n c /(n H+ + n H o) = 3 x 10~ 4 . 
The collisional de-excitation rate coefficients of C II with 



neutral hydrogen and electrons are represented by 7 
7 e , respectively. According to Wolfire et al (1995), 



7 H = 8.86 x 10" 10 cm 3 s _1 , 
7 e = 2.1 x I0 _7 (T/I00) _o ' 5 [I 



+ 0.484 

4\3i 



x (T/l 4 ) + 4.0I(T/I0 4 ) 2 - 3.39(T/10 



and 



(15) 



(16) 



In writing energy equation, we implicitly assumed that 
the equilibrium of the abundances of charged dust particles 
is established. Not only the ionization level, but the electrical 
charge of the grains and eventually, cooling function, may 
change due to the different mechanisms of interactions be- 
tween the species. However, as long as the time scales of the 
ionization, recombination and the grain charging are much 
shorter than the cooling time scale, it is safe to neglect such 
complicated processes. In fact, The basic set of equation (1)- 
(13) is structurally similar to the single fluid ideal Hall MHD 
equations. This structural similarity occurs because of the 
collision which glues the weakly ionized medium together as 
a single fluid (e.g., Pandey and Wardle 2008) and therefore, 
the direct effect of recombination (which occurs on the col- 
lisional time scale) on the cooling function can be ignored 
(see also, Draine and Sutin 1987). We are interested to study 
thermal instability in a typical dusty H I region. For such a 
system it has been shown by Stiele, Lesch & Heitsch (2006) 
and Ciolek & Mouschovias (1993) that the cooling time scale 
is larger than the other important time scales. 

Finally, we can write equation of state as 

P = ~PT, (17) 
A* 

where R is the gas constant and p represents the molecular 
weight. 

Equations (8), (9), (12), (13) and (17) along with the 
equation 



V.B = 0, 



(18) 



are our basic equations for magneto-thermal instability in a 
partially ionized medium with charged dust particles. 



2.2 Linear perturbations 

The set of the equations can be linearized in the absence 
of gradients in the initial state. We assume that the ini- 
tial equilibrium state are characterized by the values po, 
Po, To and Bo, and the velocity is assumed to be zero in 
the nonperturbed initial equilibrium state of the fluid and 
n(po,T ) = 0. 

We assume perturbations of the form 



n(r, t) = III exp(cot + ik ■ r) 



(19) 



where Hi is the amplitude of the perturbations, u> is the 
growth rate of the perturbations and k is the wavenumber 
of the perturbations. A condensation mode occurs when to 
is real and positive. If ui is real and negative, the perturba- 
tions damp. But when u> is a complex number, the system 
is oscillatory growing or damping depending on the sign of 
K(w). The system is called overstable, if we have K(u>) > 0. 




OS 0.10 - 




Figure 1. Growth rate of magneto-thermal condensation mode 
versus wavenumber of the perturbations when a = 1, 7 = 5/3, 
fT/fp = 1/2 and o p ok = 0.01. Dashed line represents growth 
rate in non-magnetized case and without charged dust particles. 
Each curve is labeled by the ratio yfX = fc p c s /w mc( j. 



We mainly interested in the possible effects of the charged 
dust particles on the condensation modes. 
Then, the linear equations are 



wpi + ipok ■ vi = 0, 

k Bi 

wpovi + ikPi + i(B • Bi)- i(k ■ B ) — 

47T 47T 



0, 



-Pi 



W7P0 



7 - 1 " ' (7 - l)po 
+{K\\k\ + K±kl)T 1 =0, 

wBi + iB (k ■ vi) - i(k ■ B„)vi 

c i + e 



pi + pofippi + pqQtTi 



4ir Zend 
Pi pi Ti 



(k-B )(kxBi) =0, 
= 0. 



Po Po Po 

Note that the derivative Q p = (dQ/dp)r and Qt 
(8Q/dT) p are evaluated for the equilibrium state. 



(20) 
(21) 

(22) 

(23) 
(24) 
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2.3 Dispersion relation 



We introduce the coordinate system e x , e y , and e z specified 
by 

Bo Bp x k 



Gy X 6 2 



(25) 



(27) 



* |B„ xk| ' 
Also, we introduce the following wavenumbers 
k P = M(7 - l)/9ofipORc s T )"\ fc T = //(7 - I^tCRcs) -1 , 
fc ^H = [M7- l)^||] _1 (-Rc s po), 

Vl = [M(7 - l)A'±]- 1 (iia.po). (26) 

Then, we can write the dispersion equation using the follow- 
ing non-dimensional quantities, 

^ _ u> _ k p _ kr k k 

kcs k k ^ k k ,^ k%^j_ 

Thus, the characteristic equation becomes 

r 7 + p 6 r 6 + p 5 r 5 + p 4 r 4 + p 3 r 3 + p 2 r 2 
+p x r + p = o, 

where the coefficients are 

Pi = q 2 ( 2 7 _1 (ct t + OK - <7 p ), 

Pi = <v, 

P2 = a\(<j T + o- K ) - a<7 _1 (0> - or - ck) 

x[2 + aA(fc/fc p ) 2 ], 

P 3 = aC[2 + Q + aA(A:/A: p ) 2 ], 

P 4 = a(ffr + <tk)[1 + C + a<A(fc/fc p ) 2 ] 

-7 _1 (a-p - or- <tk), 

P 5 = 1 + a + «C[1 + aA(fc/fc p ) 2 ], 

P 6 =cr T +(TK, (28) 

and £ = cos 2 # and # is the angle between Bo and k and 
a = (va/c b ) 2 where va is the Alfven velocity. Also, we have 
UK — C-ffyC + a K ± (l ~ C)- However, we note that thermal 
conductivity in a plasma is suppressed perpendicular to the 
magnetic field if the electron gyroradius is much smaller than 
the collisional mean free path (Spitzer 1962), as is gener- 
ally true under interstellar conditions. Also, we have A = 
(fcpCs/oJmcd) 2 , where w mc d = 0»d/po)(l/l + @)(ZeB/m d c) 
is the modified dust-cyclotron frequency and (fepC s ) _1 rep- 
resents the cooling time-scale. Thus, parameter yfX is ra- 
tio of the cooling rate and the modified dust-cyclotron fre- 
quency. We see that possible effects of the charged dust 
particles on the unstable modes of equation (27) appears 
through parameter A. If we set A = 0, our dispersion equa- 
tion (27) reduces to the standard characteristic equation of 
the magneto-thermal instability, but without dust particles 
(Field 1965). In this case, dispersion equation (27) becomes 



(r 2 + aC){r 5 + (<7 T + a- K )r 4 + (1 + a)r 3 + [a(a T 



- 0"k) 




+ (CT K + CTT-CT p )/7]r 2 +aCr + aC(CTT + CTK-ap)/7} = 0.(29) 

The expression inside the above parenthesis describes stable 
modes. But fifth order polynomial that has been obtained 



1 10 
a (Angstrom) 



Figure 2. Modified dust-cyclotron frequency u> m cd versus radius 
of the dust particles (in Angstrom) with various electrical charge. 
We assume the level of the ionization is Pi/po = 10~ 6 and the 
magnetic field strength is Bq = 10~ 6 G. Also, po/m p = 100 cm -3 
and pd/po = 0.01. Each curve is labeled by \Z\. 



by many authors (e.g., Field 1965; Stiele et al 2006) gives us 
magnetothermal condensation modes without charged dust 
particles. 



3 ANALYSIS 

Now, we can study the effects of charged dust particles 
on the condensation modes in typical plasma containing 
charged dust particles. Equation (27) is of odd degree in T 
and must therefore admit at least one positive real root for 
non-perpendicular perturbations (i.e., ( =fc 0) if the last term 
Po is negative. But when £ becomes zero (i.e, wavevector is 
transverse to the magnetic field) the dispersion equation (27) 
reduces to the classical case of the transverse perturbations 
(see Field 1965, Eq. 61). In this case, Hall like term (due to 
the existence of the charged dust particles) in equation (12) 
does not contribute to the problem and the criterion of the 
instability becomes P4 < that is consistent to the previous 
studies (e.g., Field 1965). Thus, we see that the condition 
for a positive real root (implying monotonic instability) for 
a transverse field is independent of the presence of charged 
dust particles. Note that in our analysis, we assumed that 
dust particles do not contribute to the net cooling of the 
system. The other roots may be complex conjugates of each 
other correspond to either damped waves or growing waves. 
If all roots of equation (27) have negative real parts, then 
the system is stable magneto-thermally. Using Hurwitz anal- 
ysis, we determined conditions under which all roots have 
negative real parts. We noticed that conditions of the stabil- 
ity are the same as magnetized case without dust particles. 
In other words, dust particles do not change criteria of the 
stability of the system. 

In the absence of dust particles, the maximum effect 
of magnetic field occurs for perturbations perpendicular to 
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Figure 3. Growth rates of the magneto-thermal unstable mode 
in the typical dusty H I region, where T = 100 K, foil = 0.01, 
po = 100 cm" 3 , p-Jp = 10" 6 , p d /p = 0.01 and B = 10" 6 . 
Also, we have a = 0.048, 7 = 5/3, <tt/(T p = 0.92. Black curves 
are corresponding to a case with <t p <tx = 2.19 X 10 -6 , and for 
comparison, red curves are for a p apc = 2.19 X 10 — 4 . Solid curves 
are representing unstable modes without charged dust particles. 
Radius of dust particles is a = IA and the electrical charge is 
\Z\ = 1. It is evaluated k p c B = 9.28 X 10~ 14 s and w mcd = 6.68 X 
10- 20 s. 

the initial magnetic field, i.e £ = 0. For such perturbations, 
dust particles do not change growth rate of the condensation 
modes because all terms including A vanish when ( = 0. 
So, nearly perpendicular perturbations are considered in our 
analysis. 

We solve numerically the roots of dispersion equation 
(27) for some values of the nondimensional parameters and 
taking 7 = 5/3. We take the parameters a = 1, <tt /<j p — 1/2 
and (7 p (j k = 0.01 for comparison to Field (1965). Figure 1 
shows growth rate of the condensation modes for nearly per- 
pendicular perturbations, i.e. £ = 0.1 and 0.01. Each curve is 
labeled by corresponding ratio k p c s /uj mc d = VA. For com- 
parison, growth rate of the non-magnetized condensation 
mode without dust particles is shown by dashed lines in 
Figure 1. We can see that growth rate in the magnetized 
case without dust particles (i.e. A = 0) is lower than the 
growth rate in the non-magnetized case. Charged dust par- 
ticles enhance growth rate of the condensation mode and 
parameter A has a vital role. When cooling rate is smaller 
than the modified dust-cyclotron frequency, enhancement of 
the growth rate due to the dust particles is negligible. But 
as the cooling rate becomes larger than the modified dust- 
cyclotron frequency, growth rate of the unstable mode sig- 
nificantly increases and tends to the non-magnetized growth 
rate when parameter A is large. 



Profile of the growth rate reaches to a maximum value 
for a wavenumber fc max which depends on the input param- 
eters. Wavenumber fc max corresponding to the maximum 
growth rate increases because of the magnetic field in the 
absence of the dust particles. In other words, the most un- 
stable mode occurs at a smaller wavelength because of the 
magnetic field. But wavenumber fc max decreases because of 
the charged dust particles. Note that growth rate of the per- 
turbation parallel to the initial magnetic field line does not 
modify because of the charged dust particles. But as pertur- 
bations tend to the perpendicular direction to the field line, 
the influence of charged dust particles on the growth rate 
becomes more significant. 

Figure 2 shows the modified dust-cyclotron frequency, 
w mc d, versus the grain radius a for different electrical charge. 
We assumed the number density of the neutral component 
and the magnetic field strength are about 100 cm" 3 and 
about 10~ 6 G, respectively. Also, the ionization degree is 
assumed to be Pi/po = 10~ 6 . We see that the modified dust- 
cyclotron frequency decreases with the electrical charge. 
Based on the definition of the modified dust-cyclotron fre- 
quency, we see that oj mc( j oc \Z\/Q. On the other hand, we 
have 9 oc \Z\ 2 . Thus, we can conclude t^ mc d oc \Z\ . For in- 
creasing values of | Z\ the modified dust-cyclotron frequency 
<^mcd decreases which implies less magnetic coupling and 
larger A for a fixed cooling rate. Since growth rate increases 
with A, we can conclude that a larger electrical charge \Z\ 
implies greater growth rate. We note that the dust grains 
can only respond to variations of the magnetic field with 
frequency less than u m cd- 

Dependence of the dust-cyclotron frequency on the 
grain radius is interesting. While for the grains with radii 
between 1 and 30 Angstrom the modified dust-cyclotron fre- 
quency is constant more or less, we see that aj mc d increases 
with the grain radius for a < lA, and, w mc d decreases with 
the grain radius for a > 30A. So, the level of magnetic cou- 
pling of the charged dust particles highly depends on the 
grain radius. Since w mc d increases with a for the grains with 
radii smaller than lA, the particles are more well-coupled 
magnetically and moreover parameter A decreases for a fixed 
cooling rate which implies smaller growth rate according to 
Figure 1. Minimum influence of the grains on the growth rate 
occurs for the grains with radii between 1 and 30 Angstrom, 
for which dust-cyclotron frequency reaches to a maximum 
value. In this case, growth rate is independent of the grain 
radius more or less as long as the cooling rate is kept fixed. 
But as the grain radius becomes larger than 30A, the fre- 
quency oJmcd decreases with radius a and the dust particles 
become less coupled to the field lines. Then, parameter A 
increases with the radius a and the growth rate increases 
according to Figure 1. 

Now, we can study magneto-thermal instability in a 
typical H I region with charged dust particles. Our input 
parameters are similar to Fukue & Kamaya (2007) who 
analyzed thermal instability in a typical H I by consider- 
ing ions and neutrals, separately. But they neglected dust 
particles. Our analysis is restricted to weakly ionized case 
and the bulk density and velocity are determined by the 
neutrals. Figure 3 shows the growth rates of the unstable 
modes when the ionization degree of a typical H I with 
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temperature around 100 K and the density 100 cm~ 3 is as- 
sumed to be 10 -6 . The strength of the magnetic field is 
Bo = 10~ 6 G. Having the net cooling function (14), we 
can simply calculate that k p c s = 9.28 x 10~ 14 s. Also, we 
have feu = 0.01, p d /p = 0.01, a = 0.048, 7 = 5/3, 
ct I o p — 0.92. Black curves are corresponding to a case with 
OpOK = 2.19 x 10~ 6 , and for comparison, red curves are for 
Op(JK = 2.19 x 10~ 4 . Solid curves are representing unsta- 
ble modes without charged dust particles. Radius of dust 
particles is a = 1A and the electrical charge is \Z\ = 1. 
We can also simply calculate that w mc d = 6.68 x 10~ 20 s. 
Figure 3 shows that the existence of the charged dust par- 
ticles enhances the growth rate of the unstable modes. This 
enhancement is more significant as thermal conduction be- 
comes more effective. Also, because of the charged dust par- 
ticles the maximum wavenumber fe max shifts to the smaller 
values. It implies that in a typical H I region one may expect 
formation of larger structures due to the thermal instabil- 
ity with charged dust particles. Not only the structures are 
bigger, but they are forming faster comparing to a similar 
system without charged dust particles. 

Although Fukue & Kamaya (2007) did not give a simple 
physical explanation for their findings, we think stabilization 
of ions via drag force is understandable and it will help us 
to explain our results. First, we consider a two-fluid (i.e., 
ion and neutral) system without magnetic field. Direction 
of the exerted drag force on each charged species is oppo- 
site to the direction of its velocity. So, it acts as a restoring 
force in response to any perturbations. Obviously, the effect 
is much stronger for ions because of their mean molecular 
weight comparing to the neutrals as was shown in Figure 1 
of Fukue & Kamaya (2007). This effect is survived even in 
the presence of magnetic field lines. Because the magnetic 
tension vanishes for purely transverse perturbations and the 
magnetic pressure adds up to the thermal pressure which 
implies a more thermally stable system. In our analysis, as 
we discuss, the drag force is not significant comparing to the 
magnetic force at least for ions and electrons. So, response 
of the system to the perturbations is mainly determined by 
the distribution of the magnetic field lines. There is a drift 
velocity between charged dust particles and the neutrals. 
This fact is shown by a modified induction equation with 
an extra term similar to the Hall term. Therefore, we can 
expect dissipation of the magnetic field and reduction to the 
magnetic pressure due to the presence of charge dust par- 
ticles. Charged dust particles are well tied to the magnetic 
field lines, i.e. A ~ 0. When this parameter increases due to 
the reduction of w mc d (for a fixed cooling rate), the dust par- 
ticles are less coupled to the field lines comparing to a case 
with A ~ 0. Thus, the coefficient of the Hall term increases 
that is proportional to the inverse of w mcd . 

There are also more points before comparing our re- 
sults to Fukue & Kamaya (2007). In our approach that is 
based on Pandey & Vladimirov (2007), ions and electrons 
are strongly tied to the magnetic field lines and the den- 
sity of the charged species are neglected comparing to the 
bulk density of the system that is controlled by the neutrals. 
Also, bulk velocity of the system is determined by the neu- 
tral particles. But analysis of Fukue & Kamaya (2007) is 
within ambipolar regime by considering the complete sets of 



equations for a two-fluid (i.e., ion and electron) system. In 
other words, although we are studying a multifluid system 
(including charged dust particles), our main equations are 
still single fluid equations. More importantly, our equations 
are valid as long as ions and electrons are coupled to the 
magnetic field lines. It implies that the drag force is negligi- 
ble comparing to the electromagnetic force. For this reason 
collision does not seem to play a significant role in dynam- 
ics of ions and electrons. On the other hand, behavior of 
charged dust particles to the perturbations is understood 
by the modified cyclotron frequency. There is also another 
important point: Our result does not contradict Fukue & 
Kamaya (2007) regarding to the reduction of the growth 
rate due to collision. Actually, our results are supportive 
to them. Let's explain how. First of all, Fukue & Kamaya 
(2007) mentioned (in the abstract of their paper) that the 
instability is suppressed via the friction. But it is true for 
ions as they clearly discussed in section 4.2.3 of Fukue & 
Kamaya (2007). So, significant reduction in growth rate is 
actually for ions (not neutrals). But in our analysis, we are 
employing a single fluid approach and what we observe as 
(in)stability is for neutrals. 



4 CONCLUSION 

We studied magneto-thermal instability in the presence of 
charged dust particles with little mobility. Although dust 
particles may contribute to the net cooling function, their 
contribution can be neglected for some ranges of tempera- 
ture (e.g., Draine & Salpeter 1979). The relative drift be- 
tween the plasma and the charged dust grains lead to the 
Hall term in the resulting induction equation. It is know 
that this feature is a generic property of the multifluid plas- 
mas (e.g., Wardle 1999; Pandey & Vladimirov 2007). The 
inclusion of charged dust particles does not alter the crite- 
rion for magneto-thermal stability, but the growth rate of 
the unstable modes are significantly modified. 

We showed that the growth rate of the condensation 
mode depends on the ratio of the cooling rate and the mod- 
ified dust-cyclotron frequency. Growth rate increases with 
the ratio. Thus, as charged dust particles become less well- 
coupled to the field lines, the growth rate becomes larger 
and the system is more stable to the perturbations as long 
as the cooling rate does not change. Also, the growth rate 
increases with the electrical charge. While the growth rate 
increases with the grain radius when they are larger than 30 
Angstrom, we see opposite behavior when grains are smaller 
than 1 Angstrom. Magneto-thermal condensation mode is 
independent of the grain radius, when they are between 1 
and 30 Angstrom. Note that these critical sizes may change 
for the another set of input parameters, but the profile of 
^mcd is something similar to Figure 2 for the other input 
parameters. 

Our analysis shows that even a small amount of the 
charged dust particles may modify magneto-thermal con- 
densation modes. However, a more detailed analysis is 
needed which should also take into account the dynamics 
of the dust and neutral particles. 
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